Analysis of the Q 2 dependence of charged-current quasielastic processes 

in neutrino-nucleus interactions 



O 

(N 

(N 



O 
=5 



(N 
> 

00 

q 

o 
o 



Artur M. Ankowski, 1, 2 '0 Omar Benhar, 2 and Nicola Farina 2 

1 Institute of Theoretical Physics, University of Wroclaw, 
Wroclaw, Poland, 
2 INFN and Department of Physics, "Sapienza" Universita di Roma, 
1-00185 Roma, Italy 
(Dated: July 14, 2010) 

We discuss the observed disagreement between the Q 2 distributions of neutrino-nucleus quasielas- 
tic events, measured by a number of recent experiments, and the predictions of Monte Carlo sim- 
ulations based on the relativistic Fermi gas model. The results of our analysis suggest that these 
discrepancies are likely to be ascribable to both the breakdown of the impulse approximation and 
the limitations of the Fermi gas description. Several issues related to the extraction of the Q 2 distri- 
butions from the experimental data are also discussed, and new kinematical variables, which would 
allow for an improved analysis, are proposed. 
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I. INTRODUCTION 

A number of recent experiments reported a sizable 
disagreement between the measured Q 2 distributions 
of neutrino-nucleus scattering events and the results of 
Monte Carlo (MC) simulations The largest dis- 

crepancies occur in the region of low Q 2 , typically Q 2 < 
0.2 GeV 2 , where the observed number of charged-current 
quasielastic (CCQE) events is significantly lower than 
MC predictions. To account for this differences, the Mini- 
BooNE Collaboration introduced an additional parame- 
ter in the relativistic Fermi gas (RFC) model used in the 
MC code @. This procedure leads to a modification of 
the treatment of Pauli blocking, whose physical interpre- 
tation is quite questionable. In Ref. [2j it was also argued 
that a more refined treatment of nuclear effects may be 
needed. However, the results of Ref. 4] suggest that re- 
placing the RFG model with the approach based on the 
use of a realistic spectral function (SF) does not lead to 
a consistent description of low- and high-Q 2 data. Note 
that, although strictly speaking the spectral function can 
also be defined in the RFG model, in which case it re- 
duces to a collection of (S-function peaks, in the following 
we will use the acronym SF to indicate a spectral function 
obtained from more advanced dynamical models. 

The basic assumption underlying MC simulations is 
the validity of the impulse approximation (IA), imply- 
ing that the scattering process involves only one nucleon, 
while the remaining (A — 1) particles act as spectators. 
This scheme is likely to be applicable if the space resolu- 
tion of the incoming neutrino, A ~ |q|~ , where q is the 
momentum transfer, is small compared to the average 
separation between nucleons in the target nucleus. 

In this paper we analyze possible sources of the ob- 
served discrepancies, associated with both the breakdown 



of the IA and the limitations of the RFG model, and 
point out some issues related to the extraction of the Q 2 
distributions from the experimental data. 

In Section [TT| we discuss the bias of the Q 2 reconstruc- 
tion from the kinematics of the CCQE events, while Sec- 
tion Mil is devoted to a quantitative investigation of the 
limits of applicability of the IA scheme. In Section ITVl we 
discuss the main difficulties involved in the comparison 
between theoretical calculations and experimental data, 
and put forth the proposal of new kinematical variables, 
whose use may allow for a more effective data analysis. 
Finally, in Section [Vj we summarize our results and state 
the conclusions. 



II. UNCERTAINTIES IN THE Q 2 
RECONSTRUCTION 

Consider the yield of CCQE events averaged over the 
MiniBooNE neutrino flux. We have compared the dis- 
tributions of events plotted as a function of Q 2 and its 
reconstructed value, defined as [l|, Hj 



Q 2 cc = -m 2 +2E^(E^-\k'\cost 



(1) 



where m M , k', and denote the muon mass, momentum, 
and scattering angle, respectively. In the above equation 
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with AM 2 = M 2 - M 2 , M n and M p being the neutron 
and proton mass, is the reconstructed energy of the in- 
coming neutrino. Note that Eqs. (JTJ) and have to 
be regarded as definitions of two quantities used in data 
analysis and should not be identified with the true Q 2 
and neutrino energy. 

The results of the calculations, carried out for a carbon 
target using both the RFG model, with Fermi momentum 
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FIG. 1. (color online). Distributions of CCQE events in car- 
bon as a function of Q 2 and Q 2 CC . The calculations have been 
carried out for the MiniBooNE flux, using the RFG model 
(upper curves) and the SF approach (lower curves). 



Pf = 220 MeV and separation energy e = 34 MeV, and 
the SF approach, with the spectral function of Ref. @, 
are shown in Fig. [TJ In the SF approach Pauli blocking 
has been taken into account using realistic nuclear matter 
momentum distributions and the local density approxi- 
mation, within the scheme discussed in Section HID It 
is apparent that in both cases the definition of Q 2 ec in 
terms of the measured quantities |k'| and 9 allows one to 
reproduce quite well the results obtained using the true 
Q 2 . However, as the results shown in Fig. Q] involve a 
flux average, one may still ask the question whether the 
two quantities, Q 2 and Q 2 CC , are totally equivalent. 

To clarify this point, in Fig. [2] we show a compari- 
son between the differential cross sections da/dQ 2 and 
da/dQ 2 cc obtained from the RFG model at fixed neu- 
trino energies E v = 1.2 GcV, 0.8 GeV, and 0.4 GeV. 
With the exception of the lowest energy, the maxima 
are again in very good agreement. For this reason the 
cross sections averaged over a not-too-low-energy flux, 
like MiniBooNE's, are in good agreement. 

On the other hand, the cross sections show quite a rich 
structure, exhibiting bumps, dips, and knees, which make 
the Q 2 and Q 2 ec distributions clearly different. While 
this structure is somewhat emphasized within the RFG 
model, its origin is largely model independent, as it 
can be traced back to kinematics. In this context, the 
most relevant feature is the location of the single parti- 
cle strength, expressed by the average separation energy, 
which turns out to be quite close in the RFG and SF 
models. 

In spite of the fact that the structure is not visible in 
Fig. [I] as it is washed out by the flux average, the results 
shown in Fig. [2] imply that neutrinos of fixed energy con- 
tribute in a slightly different manner to da/dQ 2 and to 
da/dQ 2 cc . At neutrino energy 0.8 GeV, corresponding to 
the peak of the MiniBooNE flux, the difference between 
the two cross sections is pronounced. 

Before discussing why the replacement of Q 2 with Q 2 CC 
leads to a significant change of the cross section, let us fo- 
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FIG. 2. (color online). Comparison of the differential cross 
sections da/dQ 2 and da/dQ 2 cc for neutrino energy 1.2 GeV 
(top panel), 0.8 GeV (middle panel), and 0.4 GeV (bottom 
panel) calculated within the RFG model. 



cus on the mechanism responsible for the observed struc- 
ture in the case of da/dQ 2 . To find the value of Q 2 cor- 
responding to the maximum we note that, due to Pauli 
blocking, the cross section may increase with momentum 
transfer up to the value |q|* = 2pp, needed to knock out 
nucleons sitting at the bottom of the Fermi sea. The 
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FIG. 3. (color online). Upper panel: curves of fixed Q 2 in 
the (co, |q|) plane, superimposed to the kinematically allowed 
region (within the dotted boundaries) for 0.8-GeV muon neu- 
trino scattering. Lower panel: same as in upper panel but for 
fixed values of Q? ec . 



energy transfer can be written as 



to — E p i — Ep 



VM 2 + (p + q) 2 - VM 2 + P 2 + e, 



(3) 



where M = (M n + M p )/2, while E p and p (E p > and 
p' = p + q) are the neutron (proton) energy and mo- 
mentum. The above equation shows that, for any given 
|p|, knocking out nucleons with momenta parallel to q 
requires the highest energy. As a consequence, the mo- 
mentum and energy transfer |q|* and w», with 



UJ, = y/NP + (3p F ) 2 -^M 2 +p% + 



(4) 



correspond to the maximum of the Q 2 distribution. The 
position of the maximum, Q 2 = 0.146 GeV 2 , is indepen- 
dent of neutrino energy, as long as E v is high enough for 
and |q|* to lie within the kinematically allowed region, 
and not too close to its boundary. The bottom panel of 
Fig. [2] shows that for E v = 0.4 GeV this condition is not 
fulfilled. 

The knee of the cross section, particularly visible for 
E v = 0.8 GeV, results from the phase space shrinkage 
above a certain Q 2 . The boundaries of the kinematically 
allowed region in the (w, |q|) plane, determined from the 



E v - |k'| < |q| < E v + |k'|, 

Ip'MpI < |q| < Ip'I + IpI, 

0<|p| <Pf, 



(5) 



where Ik'l 
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with 



h — pf < |q| < h + pf, 
Eu - I < |q| < E v + I, 



h = \J {Ep + uj) 2 — M 2 , 



(6) 



l = J(E v - U )*-m% 



(7) 



E F 
ui = lu — e. 



As illustrated in the upper panel of Fig- El where the kine- 
matically allowed region for E v = 0.8 GeV lies within the 
dotted lines corresponding to the above conditions, the 
available phase space starts to decrease above Q 2 — 0.607 
GeV 2 . As a consequence, the differential cross section 
also decreases, exhibiting a knee starting at this value 
of Q . Comparison between the results corresponding to 
E v = 0.8 GeV and 1.2 GeV shows that as the neutrino en- 
ergy gets higher this effect becomes less significant. This 
pattern can be understood considering that for higher E v 
the reduction of the phase space starts at higher |q| and 
Q 2 , where the cross section is smaller. In addition, the 
shrinkage of the phase space is less pronounced, as the 
allowed Q 2 range is much broader. 

The relation between Q 2 and Q 2 ec is by no means sim- 
ple. For example, knowing only the energy loss u> and 
the momentum transfer |q|, in MC simulation we cannot 
determine Q 2 ec ; one additional independent quantity, e.g. 
the neutrino or muon energy, is necessary. This implies 
that using Q 2 CC we are not describing the intrinsic target 
response: some information on the interaction vertex is 
also involved. To calculate Q 2 ec we may use the equality 



E^ - |k'|cos0 



Q 2 



2E V 



following from the identity 



k ■ k' 



\[{k-k') 



(8) 



(9) 



where k = (E v ,\£) and k' — (E^,^) are the neutrino 
and muon four-momenta, respectively. When we map 
Q 2 onto Q 2 ec using Eqs. JT]) and ©, the resulting val- 
ues arc shifted as listed in Table HI which explains why 
the structures in da/dQ 2 and da/Q 2 . cc appear at differ- 
ent positions. The lower panel of Fig. [3] illustrates the 
reason of the enhancement of the cross sections' knees, 
which turns them into bumps: where the phase space for 
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TABLE I. Position of the knees or bumps of the distributions 
shown in Fig. [2] 



E v (GeV) 


u (GeV) 


|q| (GeV) 


Q 2 (GeV 2 ) 


Ql c (GeV 2 ) 
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FIG. 4. (color online). Comparison of curves of fixed Q 2 
and Q 2 CC corresponding to the knee (bump) of the differential 
cross section da/dQ 2 (da/Q 2 cc ). 

da/dQ 2 shrinks, we observe an increase of the allowed 
range of Q 2 ec ; the behavior in the two variables is thus 
completely different. 

In Fig. |4] we compare the lines of constant Q 2 = 0.607 
GeV 2 and Q 2 CC = 0.306 GeV 2 corresponding to E v = 0.8 
GeV. The results clearly show that the two variables, in 
spite of leading to very similar distributions of events, are 
quite different. It clearly appears that the Q 2 and Q 2 CC 
distributions are determined by the nuclear response in 
different regions of the (w, |q|) plane, the latter being 
sensitive to significantly lower values of the momentum 
transfer. This feature implies that the Q 2 CC distribu- 
tion receives contributions from the kinematical region 
in which the validity of the IA may become questionable. 



III. BREAKDOWN OF THE IMPULSE 
APPROXIMATION 

In order to pin down the boundary of the kinemati- 
cal region in which the IA is applicable, we have studied 
the response of nuclear matter (a translationally invari- 
ant system consisting of equal number of protons and 
neutrons subject to strong interactions only) to a scalar 
probe delivering momentum q and energy to: 

S(q, u) = ^2(0\p{\n)(n\p^\0) S(lu + E - E n ). (10) 

n 

In the above equation, the operator p q = J2 P a p+ q a pi 
a p and a p being nucleon creation and annihilation op- 
erators, describes the fluctuations of the target density 



induced by the interaction with the probe. The target 
ground and final states, |0) and are eigenstates of 
the nuclear Hamiltonian H belonging to the eigenvalues 
Eq and E n , respectively. Note that scattering cross sec- 
tion of a scalar probe (e.g. thermal neutrons scattering 
off liquid helium) can be written in the simple form 

where (da/dil)o is the elementary cross section describ- 
ing scattering off individual target constituents, while the 
response defined in Eq. (|10[) is an intrinsic property, fully 
determined by internal target dynamics. 

In the case of electron- or neutrino-nucleus scattering, 
Eq. (|10j) can be readily generalized, replacing the density 
fluctuation operator with the operators describing the 
electromagnetic and weak nuclear currents. The resulting 
response tensor reads 

W v » = ^(0|jt|n)(n|J M |0) 5 4 {q + p -p n ). (12) 

n 

In the IA regime, Eq. fjlOp can be rewritten in the form 

S(q,w) = J d 3 pdEP h (p,E)P p {p + q,u-E), (13) 

where the hole and particle spectral functions, Ph and P p , 
describe the energy and momentum distributions of the 
struck particle in the initial and final states, respectively. 
The simplest implementation of the IA (usually referred 
to as plane wave impulse approximation, or PWIA) is 
based on the further assumption that final state inter- 
actions (FSI) between the knocked out nucleon and the 
spectator particles can be neglected. As a consequence, 
the nuclear matter particle spectral function reduces to 

P p {p + q,u-E) = [l--irp% n(p + q)] 

x 5(u - E + M - E^), (14) 

where n(p + q) is the occupation probability of the single 
particle state of momentum p + q and E p denotes the en- 
ergy of a free nucleon carrying momentum p. Note that, 
in the above definition, Pauli blocking of the phase space 
available to the struck nucleon is treated in a consistent 
fashion, using the momentum distribution obtained from 
the hole spectral function through 

n(p) = J dEP h {p,E). (15) 

In most calculations of electron- and neutrino-nucleus 
scattering cross sections, Airppn(p) / 3 is replaced with 
the RFG result 9{p F — |p|), yielding 

P p (p',£) = [l-&(p F -\p'\)]S(£ + M-E p ,) 
= 6(\p'\ - PF )6(£ + M-E P >) . 
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It has to be pointed out that the above prescriptions, 
while being justified in the case of uniform nuclear mat- 
ter, are questionable when applied to nuclei. In nuclear 
matter, due to translation invariance, the linear momen- 
tum is a good quantum number, that can be used to label 
single particle states. As a consequence, the momentum 
distribution also provides the occupation probability of 
the states. In finite nuclei, on the other hand, single 
particle states must be labeled according to the total an- 
gular momentum J. In this case, for any given p, n(p) 
receives contributions from states of different J, and may 
even exceed unity. 

The available results of accurate nuclear matter calcu- 
lations can be used to model the particle spectral function 
of finite nuclei within the framework of the local density 
approximation Q, i.e. using the definition of Eq. (TTJ} 
with 

^np F n(p) -> Jd 3 r^Trp 3 F n N M{p{r),p]pA(r), (17) 

where pa{t) is the nuclear density distribution, normal- 
ized to unity, and nNM[p{r), p] is the momentum distri- 
bution of nuclear matter at uniform density p(r). This 
procedure has been used in all calculations of nuclear 
cross sections discussed in this paper. 

Within the nonrelativistic approximation, in which 
both the response and the hole spectral function can be 
evaluated using realistic nuclear Hamiltonians, the valid- 
ity of the IA can be tested comparing 5(q, u) of Eq. (fTQ|) 
to 
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_600 MeV 
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5p W iA(q,w) = J d 3 pdEP h {p 1 E) 



4?r , , ' 



x5\u)-E- 



|p + q|- 

2AI 



(18) 



for different values of the momentum transfer q. 

The nuclear matter S(q,ui) at equilibrium density, 
po = 0.16 fin (corresponding to pf = 262 .4 MeV), has 
been recently computed using the correlated basis func- 
tion formalism and an effective interaction derived from 
a state-of-the-art parametrization of the nucleon-nucleon 
potential [6j. To analyze the interplay between short- 
and long-range correlations, the response defined as in 
Eq. (fl7J|) has been evaluated in both the Hartree-Fock 
and Tamm-Dancoff approximations. 

Figure [5] shows a comparison between the responses 
of Ref. [6j, obtained using Eq. (fTOj) and the correlated 
Hartree-Fock approximation, and those obtained from 
Eq. (fl"8|) using the nuclear matter hole spectral function 
of Ref. [7[ . The main difference between the two calcu- 
lations lies in the treatment of the target final state. In 
the IA scheme the state describing the struck particle is 
factored out, while in the approach of Ref. [6j the final 
A-nucleon state includes both statistical and dynamical 
correlations between the struck particle and the specta- 
tors. To make the comparison fully consistent, the PWIA 
response has been computed including only the contribu- 
tions of one-hole final states to Ph(p, E). The results of 
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FIG. 5. Energy dependence of the nuclear matter response 
S(q,uA. Solid lines: correlated Hartree-Fock approxima- 
tion [6J. Dashed lines: PWIA results, obtained from Eq. (|18p 
using the SF of Ref. 0] • Dot-dashed lines: results of the Fermi 
gas model (with nonrelativistic energies) at hp = 262.4 MeV, 
corresponding to the equilibrium density of nuclear matter. 
The panels are labeled according to the value of |q|. 



the Fermi gas model with nonrelativistic kinetic energy 
spectrum are also displayed. 

The results of Fig. [5] clearly show that at |q| < 2p F 
the response obtained from Eq. (TTSl does not exhibit 
the linear behavior at low ui resulting from the antisym- 
metrization of the final state. On the other hand, as the 
momentum transfer increases, the PWIA response draws 
closer to the one obtained in Ref. @. At |q| ~ 600 MeV 
the results of the two approaches are within 10% of one 
another in the region of the maximum. Note that inclu- 
sion of dynamical FSI, e.g. according to the approach 
of Ref. 8] , would produce a quenching of the PWIA re- 
sponse in the top panel of Fig. O thus bringing the solid 
and dashed lines in even better agreement. Theoretical 
studies of electron-nucleus scattering pl| suggest that at 
q| ~ 600 MeV, the FSI effect at the quasifree peak is 
-10 %. 

The emerging pattern suggests that the assumptions 
underlying the IA are likely to be valid at momenta larger 
than ~2p F , while at lower |q| factorization does not ap- 
pear to provide an adequate description of the final state. 
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QL (GeV 2 ) 

FIG. 6. (color online). Comparison of the MiniBooNE 
parametrization of the data (dotted line), labeled as the k- 
fit, to the spectral function calculation (dashed line). The 
solid line depicts the contribution to the latter from the re- 
gion where the IA is expected to be valid. The SF results are 
multiplied by a factor 1.12 to make them match the K-fit. 

In addition, it has to be kept in mind that at |q| ~ pp 
long-range correlations, involving more than one nucleon, 
also play a significant role @. 

Despite the fact that, being based on the nonrelativis- 
tic approximation, the approach of Ref. @ should not 
be used in the calculation of the MiniBooNE event dis- 
tribution, it helps to clarify why the Q 2 CC distribution 
of events measured by this experiment can be described 
by the RFG model [2j and the SF approach only for 
Q r 2 cc > 0.25 GeV 2 . 

To improve the description at lower Q 2 , the Mini- 
BooNE Collaboration introduced the ad hoc additional 
parameter k in the RFG model. While the authors of 
Ref. 2j argue that this procedure allows for a better 
treatment of Pauli blocking, we think that the inclusion 
of k cannot be justified on physics ground. Therefore, we 
refer to the modified RFG model of Ref. |2J as the K-fit. 

Figure[6]shows that in the region of lower values of Q 2 CC 
the contribution of |q|'s below 2pp is dominant. Hence, 
to explain the neutrino- nucleus cross section da/dQ 2 cc at 
quantitative level in the whole range of Q 2 ec one needs a 
consistent description of both short- and long-range cor- 
relations. Unfortunately, due to the difficulties involved 
in the treatment of relativistic particles in the final state, 
such a description has not been fully developed yet. 



IV. EXPERIMENTAL ISSUES 

The average axial mass extracted by past experiments 
turns out to be significantly lower than the values re- 
cently reported by K2K and MiniBooNE @. Although 
the early experiments were carried out using a variety of 
targets, about 60% of the total number of recorded events 
were collected on deuteron. The results of deuteron mea- 
surement with highest statistics [Tol - fl2T | are very consis- 
tent with one another and have small error bars. Com- 



pared to deuteron experiments, those carried out with 
heavier targets have much poorer statistics, the typical 
number of events being lower by one order of magnitude. 
Moreover, using a deuteron target allows one to mini- 
mize the systematic error associated with the treatment 
of nuclear effects. 

A lower value of the axial mass is also supported by 
the result recently obtained by the NOMAD Collabora- 
tion [l3[ using a carbon target. 

On the other hand, MiniBooNE collected more events 
than all other experiments combined, and carried out 
an analysis based on the shape of the reconstructed Q 2 
distribution. MiniBooNE reported 193709 events sur- 
viving the cuts, of which about 180000 correspond to 
< Q 2 CC < 1 GeV 2 . In the analysis not involving the 
additional parameter n, two cuts were applied: one at 
low Q 2 ec , to exclude the region where the RFG model 
is expected to break down, and one at high Q 2 oc , to ex- 
clude the region of low statistics. After rejection of the 
events excluded by the cuts, the data sample employed 
to extract the axial mass reduced to ^112000 events, i.e. 
^62% of the total. Note, however, that this figure is still 
^60 times larger than the number of events typically col- 
lected in deuteron experiments. 



A. Identification of CCQE events 

The first problem to be addressed in the comparison 
between theoretical calculations and experimental results 
is background simulation, which is, to a significant ex- 
tent, detector dependent. 

The extraction of the axial mass requires an accurate 
selection of the CCQE events and a quantitative descrip- 
tion of the irreducible backgrounds which may change 
the topology of the observed event. For example, a 
A-production process followed by pion absorption may 
be undistinguishable from a CCQE interaction, whereas 
a primary CCQE process may yield pions due to final 
state interactions. Therefore, a more complete theoreti- 
cal analysis should account for intranuclear cascade in a 
consistent way. 

Although cascade calculations including nucleon- 
nucleon correlations have not been developed yet, we can 
gauge the relevance of these effects using the results of a 
comparison of the MC generators employed in the anal- 
ysis of neutrino experiments, reported in Ref. [14]. The 
predicted ratio between true and observed CCQE events 
turns out to be in the range 82%-89% for a 1-GeV 
beam and oxygen target. On the other hand, the frac- 
tion of events misidentified as inelastic due to pion pro- 
duced in final state interactions is typically 0.5%-3%. 
The uncertainty in the modeling of pion production and 
absorption may be partly responsible for the disagree- 
ment between the values of the axial mass reported by 
K2K and MiniBooNE and those resulting from different 
experiments. 

Scattering off a correlated pair of nucleons may also 
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be misinterpreted as pion production, as it produces two 
hadron tracks in the final state. However, even neglect- 
ing nucleon absorption this background does not exceed 
3% of the CCQE events in the MiniBooNE kinematical 
range, and is therefore not likely to be significant. Be- 
cause of the weaker Q 2 dependence, resulting from the 
higher nucleon removal energy, a proper interpretation 
of these events may only marginally increase the value of 
axial mass extracted from the analysis. 

The effect of FSI, negligible in light nuclei and not 
taken into account in this paper, is expected to make the 
Q 2 CC distribution flatter, quenching its maximum and re- 
distributing strength towards higher values of the four- 
momentum transfer [l5j]. This behavior is due to the 
fact that FSI couple one-hole states to one-particle-two- 
hole states, thus leading to an increase of the average 
removal energy. As a result, the neutrino-nucleon inter- 
action takes place at higher Q 2 . Hence, part of the dis- 
crepancy between the results of K2K and MiniBooNE, 
on the one hand, and deuteron-based experiments, on 
the other hand, may be ascribable to FSI effects. 

It is also very important to realize that the influence 
of reaction mechanisms beyond the IA cannot be mini- 
mized increasing the beam energy [l6j . It turns out that 
low momentum transfers (|q| < pf) provide almost the 
same contributions to the CCQE cross section for neu- 
trino energy 0.4 GeV and 100 GeV. In the absence of a 
fully consistent theoretical description of all the relevant 
mechanisms, the most reasonable option for the experi- 
mental analysis appears to be a cutoff, to reject events 
with |q| < 2pp. 

It should be kept in mind that in the region where the 
IA is not valid the cross section may not scale with the 
number of nucleons. As a consequence, the axial mass 
(or any other parametrization of the axial form factor) 
extracted from neutrino-nucleus data at low Q 2 CC may 
turn out to be target dependent. 

B. Proposal of new kinematical variables 

The axial mass is currently extracted from experimen- 
tal data using the shape of the <2r CC distribution of CCQE 
events. However, as pointed out in Section HH Q 2 CC can- 
not be directly measured; its definition [see Eq. @] in- 
volves the nucleon separation energy and depends on the 
applied approximations, e.g. the assumption that the 
struck nucleon be at rest. 

It may be convenient replacing Q 2 CC with the new, 
model independent, variable 

= E, l -\]s!\caaO, (19) 

whose definition only involves measured quantities. From 

n k-k' Q 2 + m 2 

it follows that the j3 distribution exhibits the same be- 
havior as the Q 2 ec distribution, and can be comparably 
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P (GeV) 

FIG. 7. (color online). Distribution of the MiniBooNE data 
calculated using the "ft fit" (dashed line), the RFG model 
(dotted fine), and the SF approach (sofid fine) as a function 
of the variabfe /3, defined in Eq. (|19p . 

useful in the extraction of the axial mass or in the anal- 
ysis of nuclear effects. Figure [7] shows the distributions 
obtained for the MiniBooNE flux as a function of f3. 

An even better choice appears to be provided by the 
variable 

m M + p 

Figure [SJ showing the 4> distributions of the Mini- 
BooNE data obtained using different approaches, clearly 
illustrates that the main advantage of using cj> lies in the 
fact that the deviations of the k fit from the RFG model, 
reflecting the breakdown of the IA, show up in the high 
<j> tail of the distribution. As a consequence, a single cut 
at <fi — 3 GeV -1 allows one to reject both the region of 
low statistics and the region where effects beyond the IA 
are expected to be important. Hence, the extraction of 
the axial mass from the <j> distribution would be based on 
a larger data sample. Moreover, the dependence of the <j> 
distribution on the axial mass turns out to be only visible 
at > 1.6 GeV -1 , corresponding to the region of highest 
statistics. Note that in the case of the Q 2 CC distribution 
the situation is reversed. 



V. CONCLUSIONS 

In this paper, we have analyzed possible sources of the 
observed disagreement between the Q 2 ec distributions of 
CCQE events reported by several recent experiments and 
the prediction of Monte Carlo simulations. 

As far as the treatment of nuclear effects is concerned, 
our work suggests that, in addition to the known limita- 
tions of the RFG model, discussed in, e.g., Refs. |l^[l7|. 
the most critical feature of present analyses is the as- 
sumption that the IA scheme be applicable over the whole 
range of Q 2 . Electron scattering studies have provided 
ample evidence that the IA breaks down in the region 
of low momentum transfer, which turns out to provide a 
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<t> (1/GeV) 

FIG. 8. (color online). Same as in Fig. [7] but plotted as a 
function of the variable 0, defined in Eq. (In- 



significant fraction of the observed CCQE events, inde- 
pendent of neutrino energy. 

The failure of the IA is clearly exposed by the com- 
parison between the full nuclear matter response and the 
IA result, discussed in Section IIII1 showing that at mo- 
mentum transfers |q| < 2pp the contributions of more 
complex reaction mechanisms become important, or even 
dominant. 

While our calculations focus on the effects of factor- 
ization of the final state, it must be pointed out that dif- 
ferent mechanisms should also be considered. For exam- 
ple, meson exchange currents, which are long known to 
provide appreciable contributions to the electron-nucleus 
cross section at the quasi elastic peak and beyond, are 
also expected to contribute to the background in the case 
of neutrino-nucleus scattering. 

The development of a consistent treatment of scat- 
tering processes at low and high momentum transfer 
within a formalism easily implementable in MC simu- 
lation, while being feasible, involves severe difficulties, 



and will require a significant effort in the years to come. 
On the other hand, we believe that introducing ad hoc 
modifications of the available models, lacking a sound 
physical interpretation, will not help to clarify the origin 
of the disagreement between theoretical predictions and 
observations. 

Among the issues related to data analysis, identifica- 
tion of CCQE processes appears to be prominent. The 
discussion of Section IIV A[ based on the results of MC 
simulations, indicates that it may be at least partly re- 
sponsible for the disagreement between the values of the 
axial mass recently reported by K2K and MiniBooNE 
and those obtained from different measurements. 

Improving event identification will require a more real- 
istic description of FSI, combining the intranuclear cas- 
cade approach with a fully realistic description of the 
target nucleus, and including the relevant inelastic chan- 
nels leading to pion production. The key elements needed 
to pursue this project, i.e. in-medium nucleon and pion 
cross sections and nuclear wave functions including cor- 
relation effects, can be extracted from the available data 
and from theoretical results of accurate many-body cal- 
culations. 

Finally, we suggest that data analysis might be im- 
proved using a new kinematical variable which, unlike 
Qr ec , can be defined in terms of measured quantities only. 
In addition, using the new variable would allow one to re- 
ject both the region of low statistics and the region where 
effects beyond the IA are expected to be important with 
a single cut. As a result, the extraction of the axial mass 
would be based on higher event statistics. 

ACKNOWLEDGMENTS 

A.M. A. was supported by the Polish Ministry 
of Science and Higher Education under Grant No. 
550/MOB/2009/0. O.B. and N.F. were supported by 
INFN under Grant PI31. 



[1] R. Gran et al. (K2K Collaboration), Phys. Rev. D 74, 
052002 (2006). 

[2] A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), 

Phys. Rev. Lett. 100, 032301 (2008). 
[3] K. Hiraide et al. (SciBooNE Collaboration), Phys. Rev. 

D 78, 112004 (2008). 
[4] O. Benhar and D. Meloni, Phys. Rev. D 80, 073003 

(2009). 

[5] O. Benhar, A. Fabrocini, S. Fantoni, and I. Sick, Nucl. 

Phys. A579, 493 (1994). 
[6] O. Benhar and N. Farina, Phys. Lett. B680, 305 (2009). 
[7] O. Benhar, A. Fabrocini, and S. Fantoni, Nucl. Phys. 

A505, 267 (1989). 
[8] O. Benhar, A. Fabrocini, S. Fantoni, G.A. Miller, V.R. 

Pandharipande, and I. Sick, Phys. Rev. C 44, 2328 



(1991). 

[9] K.S. Kuzmin, V.V. Lyubushkin, and V. A. Naumov, Eur. 

Phys. J. C 54, 517 (2008). 
[10] T. Kitagaki et al, Phys. Rev. D 42, 1331 (1990). 
[11] K.L. Miller et al., Phys. Rev. D 26, 537 (1982). 
[12] N.J. Baker et al, Phys. Rev. D 23, 2499 (1981). 
[13] V.V. Lyubushkin et al. (NOMAD Collaboration), Eur. 

Phys. J. C 63, 355 (2009). 
[14] M. Antonello et al, Acta Phys. Pol. B 40, 2519 (2009). 
[15] O. Benhar, N. Farina, H. Nakamura, M. Sakuda, and 

R. Seki, Phys. Rev. D 72, 053005 (2005). 
[16] A.M. Ankowski, in Proceedings of NuFact08, 

PoS(NUFACT08)118. 
[17] O. Benhar and D. Meloni, Nucl. Phys. A789, 379 (2007). 



